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Abstract. The dynamical activity K{t) of a stochastic process is the number of times 
it changes configuration up to time t. It was recently argued that (spin) glasses are at 
a first order dynamical transition where histories of low and high activity coexist. We 
study this transition in the one-dimensional contact process by weighting its histories 
by exjp{sK{t)). We determine the phase diagram and the critical exponents of this 
model using a recently developed approach to the thermodynamics of histories that is 
based on the density matrix renormalisation group. We find that for every value of the 
infection rate, there is a phase transition at a critical value of s. Near the absorbing 
state phase transition of the contact process, the generating function of the activity 
shows a scaling behavior similar to that of the free energy in an equilibrium system 
near criticality. 

Keywords: Phase transitions into absorbing states (theory). Fluctuations (theory). 
Stochastic processes (theory), Density matrix renormalisation group calculations. 
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1. Introduction 

While phase transitions in equilibrium can only occur in two or higher dimensions, it is 
by now well established that they can be present in non equilibrium steady states even in 
one dimension [1]. Examples can be found in driven lattice gases and in models showing 
a transition out of an absorbing state. In the latter case, the system changes from an 
absorbing state to an active one as a parameter is varied. In the absorbing phase, the 
system is frozen in a particular configuration which it never leaves: its dynamical activity 
K{t), which equals the number of configuration changes up to some time t, becomes 
therefore zero (after relaxation to the absorbing state). In the active phase, the system 
changes state regularly and K{t) becomes nonzero. In this paper, we investigate the 
statistical properties of the activity in the contact process, a well known process showing 
an absorbing state phase transition. We do this within the formalism of thermodynamics 
of histories which was recently adapted to Markov processes [2]. 

Our study was motivated by recent work on glasses. Indeed, it has been argued 
that they show 'dynamical coexistence' between histories of high and of low activity 
[HI H] and that they therefore are at a dynamical first-order transition. This picture was 
confirmed in a recent study [5l E] of one dimensional kinetically constrained models of 
glasses [7]. These models have a dynamics with detailed balance and therefore show no 
equilibrium transition in one dimension. Yet, when their histories are weighted with 
gsK(t)^ clear evidence for a first-order dynamical transition at s = was found. More 
recently, a similar transition was found in more realistic, three-dimensional, atomistic 
models of glass formers [8] where it occurs at a small but nonzero value of s. In a spin- 
glass model, the transition occurs at s = below some temperature but at a non-zero 
s above that temperature [9]. Here we present evidence that in this respect, the one 
dimensional contact process behaves like a higher dimensional glass or spin glass. We 
find evidence for a dynamical transition at s > when the contact process is in the 
absorbing state, and at s = in the active phase. Near the absorbing/active transition 
the cumulant generating function of the activity show an interesting scaling behaviour. 
We obtain our results using a density matrix renormalisation group (DMRG) approach 
to the thermodynamics of histories that we introduced recently |10j . 

While the contact process is not a model for a glass, an extended version of the 
Fredrickson- Andersen model of glasses, shows a transition in the same universality class 
as the contact process [1] (the directed percolation class [H]). Our results on the scaling 
of the activity may therefore also hold in glass models. 

Besides its relevance for glasses, the dynamical activity has in recent years gained 
importance in non equilibrium statistical mechanics [2], [T2] since it appears in several 
contexts, such as in a generalization of the fluctuation-dissipation theorem to non 
equilibrium steady states [13]. The activity is symmetric under time reversal and 
therefore plays a complementary role to entropy production which is anti-symmetric 

This paper is organised as follows. In section 2 we introduce the contact process and 



Thermodynamics of histories for the one- dimensional contact process 



3 



show how the statistical properties of its activity can be determined from the largest 
eigenvalue of a generalised generator. We also show that a dynamical transition, if 
present, should occur at s > 0. In section 3, we briefly explain the use of the DMRG 
for these kind of problems. Our results on the phase diagram are discussed in section 
4. In section 5 we present results on critical exponents and scaling. Our conclusions 
are presented in section 6. Some details on the DMRG technique used are given in the 
appendix. 

2. The activity of the contact process 

In the contact process [H], each site of a lattice A can be occupied by at most one 
particle. Introduced as a model in epidemiology, a particle is associated with a sick 
person, whereas a vacancy corresponds to a healthly individual. It is thus convenient to 
use a spin language where an empty (occupied) site corresponds to a spin up (down). 
A microstate C of the system is then given by the set {erf,? G A}. In this paper, we 
restrict ourselves to one-dimensional lattices of L sites with open boundaries. 

The dynamics of the contact process is a continuous time Markov process for which 
the probability P{C, t) that the system is in a microstate C at time t evolves according 
to the master equation 

dtP{c, t) = E Hc' ^ c)p{c', t) - w{c ^ c")P(c, t)] (1) 

C'^C 

Here w{C —>■ C) is the transition rate to go from C to C . For the contact process 
two types of transition are allowed: an occupied site can become empty with rate 1, 
while an empty site becomes occupied with rate Here C, is the number of occupied 

neighbours and A, which is the only parameter in the model, is the infection rate (within 
an epidemiological context). For further reference, we introduce the inverse lifetime r(C) 
of the state C which is given by r{C) = X]cYC^(^ ^ ^')- 

The contact process is a standard model for phase transitions out of an absorbing 
state [1]. When A < Ac, the contact process reaches an absorbing state in which all sites 
are empty. For A > Ac and in an infinite system, the model goes to a stationary state 
with a finite density p of particles. It is known from extensive numerical investigations 
that this phase transition belongs to the universality class of directed percolation jllj . 
The scaling properties of various quantities near Ac are well characterised and precise 
numerical values for the critical exponents are known, especially in one dimension [1]. 
In finite systems, the stationary state will always be the absorbing one. It is therefore 
common to add as an extra process the creation of particles at the boundary sites with 
a rate u. This addition is not believed to influence the critical behavior of the model. 
In our calculations, we choose u = 1. 

In this paper, we are interested in the dynamical activity K{t) of the contact process. 
It gives the number of times that the system has changed configuration in a particular 
realisation (history) of the process, and this up to time t. Mathematically, the statistical 
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properties of the activity can be obtained from its (cumulant) generating function 

7r(s,A,L)^ limlln(e^^W) (2) 

where the average is taken over the histories of the stochastic process. We have exphcitly 
added the L-dependence because we will later on formulate a finite size scaling theory 
for TT. Equation ([2]) shows a large mathematical similarity with the definition of the free 
energy in equilibrium in which the inverse temperature jS = l/ksT is replaced by (— s) 
and the sum over configurations is replaced by a sum over histories. Hence, the name 
thermodynamics of histories. As we will show below, ln(e*'^^*^*'') is extensive in time, in 
contrast with the free energy which is extensive in volume. 

The probability distribution of K{t) at large times can be determined from 7r(s, A, L) 
through a Legendre transformation [15], while the average activity K{\, L), its variance 
Ak{\, L), and higher cumulants can be found as derivatives of 7r(s, A, L): 

K(A,L)= lim^(ir(t)) = ^(0,A,L) (3) 
Ak{X, L) = hm 1 \{K'it)) - {Kit))'] = ^(0, A, L) (4) 



t^oo t 



ds 



The variable s can be seen as a kind of fugacity that weights histories of the contact 
process according to their activity. For s ^ 0, non-active histories are favoured, whereas 
for s ^ 0, very active histories determine vr. In this paper, we are interested in the 
behaviour of the histories of the contact process as a function of A and s. While 
the parameter s has no clear "physical" interpretation, we will show that the one- 
dimensional contact process has a phase transition as a function of s, for all A- values. 
This transition is dynamic in nature, and separates regions in parameter space where 
there is no activity from regions with activity. This allows us to see the absorbing state 
phase transition in the contact process in a broader context. Moreover, as argued in 
the introduction the behaviour found here in one dimension could be a non equilibrium 
version of a phenomenon that occurs in (spin) glasses in two or more dimensions. 

For clarity, we now repeat a standard argument [2] that shows that 7r(s, A, L) equals 
the largest eigenvalue of a matrix H{s). We therefore introduce firstly the probability 
P{C, K, t) that the system is in configuration C and has activity K at time t. Using 
([1]), we immediately find that 

dtP{C,K,t)= J2 [w{C' ^ C)P{C',K -l,t) ~w{C ^ C')P{C,K,t)] 

C'^C 

(5) 

Consequently, the discrete Laplace transform P{C,s,t) = J2k=o^'^^ ^i^^ ^y^) evolves 
according to 

dtP{C,s,t)= ^ [w{C' ^ C)e'P{C',s,t) -w{C ^ C')P{C,s,t)] (6) 

To continue, it is convenient to introduce a matrix notation as common in the so called 
"quantum" approach to stochastic particle systems [16]. We therefore introduce a set 



{e'sf - rii) + ^{ni^i + n,+i)(e's- 



Thermodynamics of histories for the one- dimensional contact process 5 

of basis vectors \C) each corresponding to a microstate C and a vector \P{s,t)) with 
components P{C,s,t) = {C\P{s,t)) . Using this notation, the set of equations can 
be rewritten as 

dt\P) = His)\P) (7) 

The diagonal elements of the matrix H{s) are equal to minus the inverse lifetimes of 
the states, while the off-diagonal elements are given by the transition rates multiplied 
by e^. For s = 0, ([7]) reduces to the master equation ([1]) and H{0) corresponds to the 
generator of the stochastic process. We will therefore refer to H{s) as the generalised 
generator. 

Using the "quantum" notation of ^16j one can easily show that for the contact 
process 

m = E 

i=l 

+ uj{s^ + sf^ -vi- vl) (8) 

(no = riL+i = 0). Here ni,Vi,sf and s~ are standard particle number, vacancy number, 
particle annihilation and creation operators at site i 

The formal solution to (j7]) is 

|P(s,t)) = e^W*|P(s,0)) (10) 
Therefore we have 

^e^m^=J2Y.^'''PiC,K,t) 

C K 

= Y.Hc,s,t) 

c 

= E(C|e^(^)*|P(^,0)) (11) 
c 

Using the spectral theorem the sum in (ITTl) can be written as a sum over the eigenvalues 
of H{s). In the long time limit in which we are interested, this sum is dominated by the 
largest eigenvalue /io(s,A,L) of the generalised generator, so that for large t we write 
am as 



= (^E(CK'o)(t'o|P(3,0))je'^»(^'^'^)Ml + ---] (12) 

Here \vq) is the eigenvector associated with /iq. Taking t — > oo the terms coming from 
smaller eigenvalues can be neglected and using ([2]) we finally obtain 

7r(s, A,L) = /io(s, A,L) (13) 

The cumulant generating function of the activity therefore equals the largest eigenvalue 
of H{s). In this way, we have also shown that ln(e*^*-*'') is extensive in time. Hence, the 
definition (121) makes sense. 
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Some insight on the behaviour of vr for negative s can be obtained from simple 
arguments. First we observe that vr is non decreasing and by definition is zero at s = 0. 
Hence we obtain the bounds 

7r(s -oo, A, L) < 7r(s, A, L) < s < (14) 

For s — oo, only histories with no activity contribute to vr. This implies that the 
system at time t still is in the microstate C in which it was initially (at time t = 0). 
The probability for this decays exponentially with waiting time l/r{C) so that 

^lim^(e^^) = ^ exp [-r{C)t] (15) 

In order to determine tt we have to take the long time limit of this expression, which 
will be dominated by min,^ r(C). Provided limits can be interchanged, we obtain 

lim 7r(s, A, L) = — minr(C) (16) 

s — >— oo C 

This result is quite general and holds also for other models. In the particular case of 
the contact process one has for the completely empty configuration, r(C) = 2uj = 2. 
For configurations with one particle r{C) = 3 + A/2 if the particle is in the bulk and 
r{C) = 2 + A/2 if it is on the boundary. Continuing in this way, one can easily see that 
mincr{C) = 2. Putting everything together we find 

-2 < 7r(s,A,L) < s<0 (17) 

This also implies that the intensive quantities tt{s, X, L)/ L and K{X,L) are constant 
and equal to zero for any negative s. Hence, if there is some non analytical behaviour 
in these quantities, it should occur at Sc > 0. 



3. DMRG approach to activity fluctuations 

To calculate the cumulant generating function of the activity of the contact process we 
have to calculate the largest eigenvalue of the generalised generator (IHl). If we interpret 
this generator as a 'Hamiltonian' [16], calculating its largest eigenvalue is mathematically 
similar to determining the ground state energy of a quantum spin chain. The main 
difference with a standard quantum problem lies in the fact that in the present case the 
'Hamiltonian' is non-Hermitian. 

One of the most succesfuU numerical procedures to calculate ground state properties 
of quantum chains is the density matrix renormalisation group (DMRG) introduced by 
S White [171 HE] (for a review see [IE]). Later, the technique was adapted for generators 
of stochastic processes [201 EI], and recently we applied it for the first time to the 
generalised generators associated with activity fiuctuations in the contact process and 
with current fiuctuations in driven lattice gases |10j . 

The application of the DMRG to H{s) is not fundamentally different from the 
standard approach used for quantum systems. For stochastic systems two modifications 
are necessary. The first one is straightforward: one needs precise diagonalisation 
algorithms capable of handling non-Hermitian matrices. These exist (e.g. the Arnoldi 
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Figure 1. Cumulant generating fmiction for A = 4.3, s < for i = 18, 20, 50. On 
the scale of the left figure, one cannot see any L-dependence. On the right, the region 
close to s = is enlarged. System size increases from top to bottom. 

algoritm [22j), but are much more time consuming than their Hermitian counterparts. 
They are also numerically less stable. This limits the sizes we can reach to L ^ 60. 

The second problem concerns the fact that left and right eigenvectors of H{s) are 
not related by transposition. In Appendix A, we briefly discuss the procedure that we 
followed to solve this problem. 

Once the largest eigenvalue has been calculated with sufficient numerical accuracy, 
the average activity and its variance are determined using numerical differentiation. 
Because of numerical errors, it is not possible to obtain cumulants beyond the third 
with high enough precision. 

4. The phase diagram 

In this section, we present evidence that there exist a phase transition between 
dynamically active and inactive phases for every A. 

Firstly, we show in Figure 1 (left side) a typical result for 7r(s, A, L). In this case 
A = 4.3, which is in the active phase (Ac — 3.2978 in one dimension). As can be seen, 
71 tends to —2 for decreasing s- values, as predicted in section 2. There is only a weak 
L-dependence, which is most pronounced near s = as can be seen on the right side of 
Figure 1. 

From results like these we can calculate the s-weighted average activity and its 
variance as a function of s for various L-values. In Figure 2, we show our results, again 
for an infection rate above the critical one. As can be seen, a peak occurs in the variance. 
The height of the peak grows with L, which is a strong indication of the occurence of a 
phase transition. The average activity has a steep increase at this point. 

The location of the maxima in the variance leads to finite size estimates Sc{L) of 
the location of the transition. Our results for these are shown in Figure 3. As proven 
in section 2, any phase transition in the infinite system must occur at some Sc > 0. 
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Figure 2. (Color online) Average s-weighted activity (left) and variance of the activity 
(right) as a function of s for different L- values (A = 4.3). The red circles indicate the 
s-values were the variance reaches its peak. 

Because of the considerable curvature that is still present when plotting Sc versus 1/L 
it is not possible to obtain very precise estimates of the location of Sc- The available 
data do indicate however, that for A > Ac, Sc = whereas for A < Ac, Sc has a small 
positive value. For example, for A = 2.3, we estimate Sc = 0.005 ± 0.001. 

In the next section, we will show that the average activity is proportional to the 
density p of particles. The (ordinary) contact process, i.e. at s = 0, always reaches the 
absorbing state (p = 0) when A < Ac- Hence, it is not surprising that one needs a finite 
value of s for the model to turn to the active state. On the other hand, when A > Ac, 
the average density is nonzero at s = 0. Since the transition has to occur at s > 0, it 
must be located at Sc = exactly. 

Other evidence for a non-zero Sc below Ac comes from a real space renormalisation 
group (RG) calculation. Since the (generalised) generators of stochastic processes can 
be interpreted as Hamiltonians of quantum spin chains, it is natural to extend real 
space renormalisation schemes introduced in the study of quantum systems [23l [211 125] 
to stochastic systems. Such an extension is possible [26] and can give very precise 
values for the critical properties of the contact process [27]. Here, we present results of 
applying this scheme also to the generalised generator (IHl). Details of this calculation, 
and applications to another absorbing state phase transition, will be presented elsewhere 
[28] . In the renormalisation we used blocks of two sites. While such a small block cannot 
give accurate numerical results for the location of critical points or for the value of critical 
exponents, it can give a good qualitative description of an RG-fiow. Our results for this 
flow are given in Figure 4. The line s = is invariant under the renormalisation. Here 
we recover the results of [26]: for A below Ac the flow is towards A = 0, whereas in 
the opposite case, A flows to infinity. The behavior becomes more interesting in the 
extended parameter space (A, s). There now exists a critical line (full red line). Starting 
from below this line, the RG flows to A = 0, whereas for parameter values above this 
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Figure 4. (Color online) Renormalisation flow for the generalised generator ([8]). The 
full (red) line separates regions that renormalise to A = from those that renormalise 
to A = oo. The other lines show renormalisation flows for various points starting near 
A — Ac, s = 0. 



line, the flow goes to infinity. It can be clearly seen, that the critical line is at s = 
for A > Ac, whereas it is at s > for A < Ac- This renormalisation calculation therefore 
supports the conclusions coming from the DMRG. 
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5. Exponents and scaling behaviour 

In the previous section we have shown numerically that for every A, the contact process 
undergoes a transition to an active state when s exceeds a threshold Sc- We now 
investigate first the behaviour of the activity at s = as a function of A. Then we 
put forward a scaling ansatz for 71(5, A, L) near the critical point at A = Ac, s = 0. 

We begin by deriving a relation between the average activity and the density of 
particles. Consider a particular history which has activity K{t) at time t. In the time 
interval between t and t + dt, the activity changes to K(t) + 1 if the model changes 
configuration which happens with probability 



P{t)dt 



L ^ L 

U=l ^ 1=1 



dt 



(here again no = '"-l+i = and we didn't take into account the effect of u). The activity 
stays constant in this time interval if the configuration doesn't change, which happens 
with probability 1 — P{t)dt. Consequently, 

d{K{t)) {K{t + dt)) - {K{t)) 



dt dt 

L ^ L 



9 -ni){ni_i +ni+i) 
u=i i=l 



The two-point correlation functions appearing in (ITS]) can be eliminated by comparing 
this equation with the equation of motion for the average particle density at site i which 
can be obtained from the master equation using standard techniques [H]. One finds 

^ = {[H{s = 0),n,]) = -{n,) + ^((1 - n,)(n,_i + n,+x)) (19) 
Eliminating the two-point correlation functions between ( |T8i) and ( |T9l) we obtain 

«.,t(n.> + i:^ (20) 

ai -^-^ -^-^ ai 

We are interested in the long time limit in which the process reaches its steady state. 
The first term on the right hand side then equals 2pL whereas the second term becomes 
zero. Hence we immediately find from the definition ([3]) 

K{X,L) = 2pL (21) 

The average activity is twice the average number of particles. 

The particle density in the contact process has well known scaling properties near 
the absorbing state transition [H [11]. Putting AA = A — Ac, one has 

p(AA, L) ~ L~f^/^^G{L^/''^AX) (22) 

where G is a scaling function and f3 and i>± are known critical exponents of the contact 
process. In one dimension, their values are f3 = 0.27649(4) and = 1.09684(6) [29i| . 
This scaling and (I2T1) imply that at Ac, the average activity behaves as a power law 

ir(Ac,L)~L'^ (23) 
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Figure 5. (Color online) Finite size estimates of the exponent a versus 1/L (for L up 
to 60). The red circle is the prediction cr = 1 — l3/vi_ = 0.7479. 




Figure 6. Variance of the activity as a function of A A for (top to bottom) 
L = 44,38,32,26,20 and 14. 

with cr = 1 — = 0.7479. We calculated finite size estimates of (t(L) = 

log[K{Xc,L + 2)/K{X^,L)]/log{{L + 2)/L) using our DMRG technique. The results 
are shown in Figure 5. Using standard extrapolation technique for finite size systems 
[30] we obtain the estimate cr = 0.746(2) in perfect agreement with the prediction made 
above. 

We now turn to the variance of the activity at s = 0. DMRG results for this 
quantity as a function of A are shown in Figure 6. The data show a clear peak whose 
height increases as a power of L. We determined the associated exponent to be 3.08(2). 

Given the mathematical similarity between the generating function vr and the 
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Figure 7. Logaritmic derivative of tlie variance (at Sc{L)) as a function of 1/L for 
various A-values. 



equilibrium free energy, it is natural to put forward a scaling ansatz for vr near the 
absorbing state transition. We propose therefore the following finite size scaling ansatz 
for 71 

7r{s, AX, L) ~ L-'F{Ly^s, L^/^^AA) (24) 

Here F{x, y) is a scaling function, and Us a new critical exponent associated with activity 
fluctuations. The exponent z is the dynamical exponent z = It replaces the 

dimension of space d that appears in the scaling of the free energy density, because vr is 
a quantity per unit of time. For the one-dimensional contact process z = 1.5805. From 
the ansatz fl24|) . we predict that the average activity at the transition scales as L~^"'"^^ 
Comparing with ( l23i) we obtain 

ys = z + l- (3/ui_ = 2.3284 (25) 

The scaling fl2^ then predicts that the variance of the activity at criticality should scale 
as L"^"*"^^" where the exponent should equal 3.0763. This value that is fully consistent 
with the numerical value found above. The DMRG results therefore fully support the 
ansatz ( !24l) . From this ansatz, one can derive the scaling behaviour of various quantities. 
For example, the quantity K{s,L) = dn/ds^s, Xc, L) should scale as L^^'^^^G^L^^'s). 
Scalings hke this are well satisfied (see figure 5 in [TO]). 

We have further investigated the behaviour of the variance of the activity along the 
transition line Sc{X). For this quantity we find a distinct behaviour below and above Ac- 
Herefore, we go back to our data on the variance of the activity as a function of s (see 
Figure 2). In Figure 7, we plot our results for d\og{d'^'7i/ds'^{sc{L)))/d\ogL for various 
L-values. The derivative is a discrete one, as obtained from a comparison of results for 
system size L and L + 2. When the variance grows as a power law, this quantity should 
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approach a constant for large L. This seems to be indeed the case for A < Ac. At Ac, an 
extrapolation leads to an exponent value of 3.05(2), consistent with the value obtained 
above from the behavior as a function of A at s = 0. An extrapolation of the data for 
A = 3 gives an exponent 2.00(2), which would imply that the variance of the activity 
per site is constant. For A = 2.3, the exponent is somewhat higher (~ 2.4). This could 
be a genuine effect, or the result of the DMRG being less stable farther away from Ac. 
In contrast, in the active region (and especially for A = 4.3), it is clear that the behavior 
of the variance of the activity is not consistent with a power law. The data can be 
fitted better with an exponential growth of the variance. Such a behavior is somewhat 
unexpected but not impossible in principle. 

6. Conclusions 

In this paper we studied the contact process and weighted its histories according to their 
activity. We studied the phase diagram of the model as a function of the parameters s 
and A and found that when s is increased the model undergoes a phase transition between 
an absorbing state with zero activity and an active state. The well known absorbing 
state phase transition at s = is thus seen in a broader context. The transition is at 
Sc > for A < Ac while it occurs at Sc = when A > Ac. Our results were obtained 
using a density matrix renormalisation group approach that allows us to calculate the 
full cumulant generating function of the activity. 

This generating function can also be obtained from simulations based on a cloning 
algorithm [311 [32]. The DMRG approach has the advantage that it gives in principle 
numerically exact results. The cloning algorithm is less accurate, but can easily be 
extended to higher dimensions, which for the DMRG is more difficult. Some preliminary 
results for the contact process using the cloning algorithm were reported in [32]. These 
authors study a somewhat different model in which particles can be created at every site 
(and not just at the boundaries as in our model). Evidence for a transition at an Sc 7^ 
was given but no systematic study of the phase diagram and the critical exponents was 
presented. 

In section 2, we showed that the s- weighted average activity (per site) is zero for 
s < Sc- Since for A > Ac, Sc = and since at s = the average activity (per site) equals 
2p, we conclude that in this range of A values the transition between non-active and 
active phases has to be first order. We are unable to give definite claims on the order 
of the dynamical transition for A < Ac. The order of the transition can be determined 
from the probability distribution of the activity. Unfortunately, it is not possible to get 
this function unambiguously from the cumulant generating function when a first order 
transition is present [33j. 

We have shown that the cumulant generating function vr has a scaling form close to 
A = Ac, s = 0. The scaling form is a natural extension of that obeyed by the free energy 
of equilibrium systems near a critical point. A similar scaling was found to hold for the 
cumulant generating function of the current in the totally asymmetric exclusion process 
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[TU] . It would be of interest to investigate whether the same scahng form apphes to the 
dynamical transition in (spin) glass models. In one dimension, kinetically constrained 
models do not show a transition (unless at T = 0) but transitions in the equilibrium 
state of glass models can occur in higher dimensions. The validity of the scaling ansatz 
in these cases would enforce the link between the contact process and models for glassy 
dynamics. 
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Appendix A. Some details on the DMRG approach 

Assume one wants to know a particular eigenstate {ip) of the Hamiltonian for a quantum 
spin chain, together with its associated eigenvalue. Often this is the ground state. 
This state of interest will be called the target state. Consider a system of size L for 
which the target state has to be determined. For a spin 1/2-system, this involves the 
diagonalisation of a matrix in a 2^ dimensional space. When L is too big this cannot 
be done exactly. The idea is then to 'project' the Hamiltonian into a space of lower 
dimension. The basis vectors spanning this lower dimensional space are determined from 
a reduced density matrix of the target vector. We now briefly describe the procedure 
used to determine these basis vectors in the case of quantum systems since this is 
necessary to explain the modifications that have to be made for stochastic systems. 

As a first step in the DMRG, the system is divided into two parts. One is called the 
system block, the other the environment. The whole system is called the superblock. 
Both parts are taken of equal size and the number of degrees of freedom for each is 
n = 2^/^. The orthonormal basis vectors representing the respective configuration space 
will be denoted by \i) and Our aim is to reduce the number of states n used to 
describe the system block without changing the target state \ip) of the superblock. More 
precisely, assume one wants to reduce n to m < n. The n vectors \i) are to be replaced 
by m vectors a = 1, . . . ,m. In combination with the environment states these 
vectors should be able to make an accurate representation {ip) of the target state {ip) 

m n 

l^)^lv;) = EE«",,>")|j) (Ai) 

a=l j=l 

By this we mean that one needs to minimise 

S = - (A.2) 

with respect to all aa,j and \u'^)- The solution of this minimisation problem is the heart 
of the DMRG and is given in terms of the density matrix p = of the target state 

of the superblock. From this matrix one constructs the reduced density matrix on the 
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system block 

n 

P^^' = Y.i^\{J\p\^')\J) (A.3) 
i=i 

p is an Hermitian operator with n eigenvalues. It can then be shown that the m 
eigenvectors of p with largest eigenvalue are the vectors \u°') in the solution of the 
minimisation problem. This result forms the basis of both the infinite size and finite 
size algorithms used in practical DMRG calculations [18]. 

When H is the (generalised) generator of a stochastic process, left and right 
eigenvectors are different. In principle there are different ways to adapt the DMRG 
approach. We followed [21] in using both the left and right eigenvector associated with 
the largest eigenvalue of H{s) as target state. In particular, if lipo) is the right eigenvector 
(and {ipol its transpose) and {ipo\ the left eigenvector (with transpose |<^o)) we used as 
density matrix the combination 

P= ^[l^o)(V^o| + l¥'o)(¥'o|] (A.4) 

This matrix is symmetric. Its reduced density matrix p is calculated as in (lA.Sp and 
the rest of the DMRG algorithms are not modified. 
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